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Precise measurement of Higgs boson couplings is an important task for International Linear 
Collider (ILC) experiments and will facilitate the understanding of the particle mass gener- 
ation mechanism. In this study, the measurement accuracies of the Higgs boson branching 
fractions to the b and c quarks and gluons, ABr(H — > bb, cc, gg)/Br, were evaluated with 
the full International Large Detector model (ILD_00) for the Higgs mass of 120 GeV at the 
center-of-mass (CM) energies of 250 and 350 GeV using neutrino, hadronic and leptonic 
channels and assuming an integrated luminosity of 250 fb _1 , and an electron (positron) 
beam polarization of —80% (+30%). We obtained the following measurement accuracies of 
the Higgs cross section times branching fraction (A(<r • Br) /a ■ Br) for decay of the Higgs 
into bb, cc, and gg; as 1.0%, 6.9%, and 8.5% at a CM energy of 250 GeV and 1.0%, 6.2%, and 
7.3% at 350 GeV, respectively. After the measurement accuracy of the cross section (Aer/er) 
was corrected using the results of studies at 250 GeV and their extrapolation to 350 GeV, 
the derived measurement accuracies of the branching fractions (ABr/Br) to bb, cc, and gg 
were 2.7%, 7.3%, and 8.9% at a CM energy of 250 GeV and 3.6%, 7.2%, and 8.1% at 350 
GeV, respectively. 

PACS numbers: 13.66. Fg, 29.20.Ej 
Keywords: ILC, Higgs boson, Branching ratio 

I. INTRODUCTION 

Precise measurement of the Higgs boson branching ratios (BRs) is an important task for the 
International Linear Collider (ILC) program. It is also crucial for the understanding of the nature of 
electro-weak symmetry breaking and provides a window to investigate physics beyond the standard 
model (SM). The relatively low background and well-defined initial state of the ILC experiments 
allow precise, model-independent study of the Higgs boson, which is not an easy task for Large 
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Hadron Collider experiments [1, 2]. Measurements of the Higgs BRs to bb and cc decays at an 
e + e~ linear collider were reported in Refs. [3-7]. In this study, we investigate the accuracies of 
BRs of the Higgs to bb, cc, and gg using Geant4 [8] based realistic simulation implemented with a 
proposed International Large Detector (ILD) [9]. 

In this study, we assume a Higgs mass of 120 GeV/c 2 and an integrated luminosity of 250 fb _1 , 
and estimate the accuracies of the BRs at center-of-mass (CM) energies of 250 and 350 GeV. The 
former value is close to the threshold of Higgs production and thus is considered as initial target 
of ILC experiments. The latter is close to the threshold of top quark pair production; therefore, 
Higgs data can be corrected simultaneously with a top threshold study. The difference between 
kinematical conditions at 250 and 350 GeV could yield different detection efficiencies and thus 
different BR accuracies. The accuracies at 250 and 350 GeV under the same conditions are studied 
and compared. 

The experimental conditions for this study are described in section II. We selected the Higgs 
events in three channels: neutrino, hadronic and leptonic. The event selection and background 
suppression processes are described in the section III. The derivation of the BRs is presented in 
the section IV, and the conclusion is given in the last section. 

II. HIGGS PHYSICS IN THE ILC EXPERIMENT 
A. ILC experiment and Higgs production 

The ILC is a future electron-positron (e~e + ) linear collider for experiments at an initial center- 
of-mass (CM) energy (\fs) up to 500 GeV, which can be extended to 1 TeV. The production cross 
section of the Higgs boson is shown in Fig. 1(a) as a function of the CM energy for a Higgs mass of 
120 GeV. At a low CM energy, the Higgs boson is produced primarily through the Higgs-strahlung 
e + e~ — > ZH process, which has a maximum around 250 GeV when the effect of the initial state 
radiation is considered. This is about 20 GeV higher than that without the initial state radiation. 
At the = 350 GeV, the total cross section is reduced, although the contribution of W/Z fusion 
is greater than that at 250 GeV. The decay BRs of the Higgs boson in the SM are shown as a 
function of its mass in Fig. 1(b). The Higgs decays mainly to bb if its mass is below 140 GeV and 
to WW* in the case of a mass of above 140 GeV. 

Higgs analysis modes are categorized in terms of the three Z boson decay channels: Z — > vv 
(neutrino), qq (hadronic), and £ + £~ (leptonic), as shown in Fig. 2. We assumed the —80% and 
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FIG. 1: (a) Production cross section of Higgs boson at a Higgs mass of 120 GeV through the Higgs- 
strahlung (ZH) (solid) and all ffH (dashed) processes assuming the -80% electron and +30% positron 
beam polarization. The cross section is calculated including the initial state radiation by Whizard [10]. (b) 
SM Higgs branching fractions as a function of Higgs mass with PYTHIA [11]. 



+30% polarization of the initial electrons and positrons, respectively, in order to enhance the Higgs 
signals. 







(a) (b) ( c ) 

FIG. 2: Higgs boson production diagrams categorized according to the final states: (a) neutrino (vvH), 
(b) hadronic (qqH), and (c) leptonic (l + t~H) channels. Each channel is produced mainly through the 
Higgs-strahlung (ZH) process at low CM energies, although the neutrino and leptonic channels also include 
the WW and ZZ fusion processes, respectively. 



4 



B. ILD concept 

We used the ILD [9] model for this study. The ILD, which is the validated detector concept for 
the ILC, is equipped with a highly segmented calorimeter and a hybrid tracking system consisting of 
gaseous, silicon-strip, and silicon-pixel trackers. They provide an excellent jet energy resolution by 
particle flow analysis, as well as excellent momentum resolution and vertex flavor tagging capability, 
which are necessary for measuring multi-jet final states in the ILC energy region. All sub-detector 
components of the ILD are shown in Fig. 3; which consists of silicon-pixel vertex detectors (VTX), 
silicon inner and outer detectors (SIT, SET), a time projection chamber (TPC), high-granularity 
electromagnetic and hadron calorimeters (ECAL, HCAL), a super-conducting solenoid magnet 
with a 3.5 T magnetic field, and an iron return yoke with a muon detector. In addition, forward 
silicon trackers (FTD, ETD) and beam/luminosity calorimeters (LCAL, LHCAL and BCAL) are 
installed in the forward region. 

6990 -i ■ 1 



Coil + Cryostat 



2020 -j= 
180£ 



bCAL 



III' 



2348 2622 3922 4072 



FIG. 3: Schematic view of the ILD detector. 



The VTX system consists of three double layers of silicon pixel sensors with a 2.8 /zm point 
resolution located at radii between 16 mm and 60 mm, the total radiation length being 0.74%. 
The impact parameter resolution (crjp) of the VTX system is 5 /im© 10 fim ■ GeV/c/psin 3 / 2 6. 
The TPC occupies a volume up to a radius of 1.8 m and a half-length in Z of 2.3 m, providing 
a stand-alone momentum resolution of cr\/p T ~ 9 x 10~ 5 GeV -1 . The SIT and SET are placed 
at the inner and outer sides of the TPC with 7 and 50 [im point resolutions in the R — <f) and 
z directions, respectively. The overall momentum resolution of the tracking system (a 1 /p T ) is 



5 



2 x 1(T 5 GeV" 1 ©1 x l(r 3 /Prsin0 for the momentum range 1-200 GeV [9]. The ECAL consists of 
24 Xq tungsten absorbers with highly segmented (5x5 mm 2 ) readouts. The HCAL consists of 5.5 
Xi steel absorbers with a 3 x 3 mm 2 scintillator tile readout. With the ILD particle flow algorithm 
package, PandoraPFA [12], a dijet energy resolution of 25% / \J E (GeV) has been achieved for a 
45-GeV dijet, which corresponds to a single-jet energy resolution of oeJEj = 3.7% [9]. 

C. Analysis framework and Monte Carlo samples 

Monte Carlo (MC) generator samples for the physics study were produced using the 
Whizard [10], and fragmentation and hadronization processes were simulated by PYTHIA [11]. The 
SM Higgs branching fractions in PYTHIA are 65.7%, 3.6%, and 5.5% for bb, cc, and gg, respectively. 
The generated particles were passed through the Geant4 [8] based detector simulator Mokka [13] 
with the ILD model. The simulated hits were digitized and then reconstructed by the MarlinReco 
package; then, the resulting skimmed data were analyzed. The statistics of the simulated Higgs 
signal samples were 500 fb _1 for both CM energies of both 250 and 350 GeV, whereas those for 
background processes varied with the signal-to-noise ratio (S/N). They are scaled in the analysis in 
order to obtain results corresponding to an integrated luminosity of 250 fb" 1 . The major SM back- 
ground processes for the e + e" — > ZH analysis are e + e~ — > ZZ and W + W~; thus we considered 
final sample states of vvqq, v£qq, Uqq, vvM, qqqq and till. In addition, the qq and ti backgrounds 
were also considered for the neutrino and hadronic channels (but only for y/s = 350 GeV because we 
used a top mass of 174.9 GeV/c 2 ). In the leptonic channel, most of the multi-jet backgrounds are 
well suppressed if dilepton identification is required; thus, only the iiqq and viqq backgrounds were 
considered. We used the 250-GeV samples produced for the ILD letter of intent (LOI) studies [9]; 
thus, their beam parameters correspond to those defined in the ILC Reference Design Report [3]. 
On the other hand, the 350-GeV samples were newly produced for this study using the updated 
beam parameter SB2009 [14]. The instantaneous luminosities were 0.75 and 1 x(10 34 cm~ 2 s -1 ) for 
250 and 350 GeV, which yield integrated luminosities of 188 and 250 fb -1 , respectively, for about 

3 years at 100 days of operation per year. 

III. EVENT RECONSTRUCTION AND BACKGROUND SUPPRESSION 

Depending on the Z decay mode, the analysis channels are categorized as the neutrino (dijet), 
hadronic (four-jets) and leptonic channels (dileptons + dijets), which are described in the following 



6 



subsections. 



A. Neutrino channel (vvH) 



For neutrino channel analysis, particles in the event are first forcibly clustered into two jets by 
the Durham jet-finding algorithm. After the dijet clustering, background reductions are applied 
according to the selection criteria in Table I. At a CM energy of 250 GeV, the Higgs is produced 
almost at rest because it is close to the production threshold, whereas it is boosted at 350 GeV. 
Thus, the cut conditions are optimized to obtain the best S/N at each energy. In this channel, Z 
boson decays invisibly (vv); thus, the vvqq and viqq processes in the SM are the main backgrounds. 
To reduce them, a cut on the missing mass (M m i ss ) is applied. Although this cut decreases the 
Higgs signal from the WW fusion process, the viqq, iiqq and qqqq backgrounds are effectively 
reduced, qq background is reduced by the following kinematical cuts: the transverse momentum 
(Pt), longitudinal momentum (Pi), and maximum momentum (P max )- The Hit background is 
well reduced by a cut on the number of charged tracks in an event (N c hd)- in addition, the viqq 
background reduction is improved by the Yyi and Y23 cuts. Y\i and Y23 are the maximum and the 
minimum of y values (scaled jet masses), respectively, required to cluster the event into two jets. 

The background reductions for each cut are summarized in Table I for each CM energy. After 
all selection criteria are met, an additional likelihood ratio (LR) cut is applied to improve the 
background reduction. The LR is defined using the following variables: M m i ss , the number of 
particles (Nppo), Y\2, Pmax, Pt, and Mjj. The likelihood cut positions are optimized to maximize 
signal significance and LR > 0.165 and LR > 0.395 are selected for CM energies of 250 and 350 
GeV, respectively. The signal significance (S/y/S + B) after all background reductions is also listed 
in Table with its efficiency, where S and B are the numbers of Higgs signal and background entries, 
respectively, after all cuts are applied. The remaining backgrounds are viqq (60%), vvqq (20%), 
and qq (10%) at both 250 and 350 GeV. 



For hadronic channel analysis, particles in the event are first forcibly clustered into four jets. 
Next, a Higgs and Z candidate dijet pair that minimize the following x 2 formula are selected from 
the four jets: 



B. Hadronic channel (qqH) 




2 



2 



(1) 
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TABLE I: Summary of the vvH channel background reduction assuming C — 250 fb 1 with P(e ,e + ) = 
(-0.8, +0.3). 



UM energy (LrCV ) 


250 


350 


Cut names 


condition 


Sig. 


Bkg. 


condition 


Sig. 


Bkg. 


Generated 




19360 


44827100 




26307 


20855900 


Missing mass (GeV) 


80 < M mlss < 140 


15466 


6214050 


50 < M miss < 240 


23202 


5627040 


Transverse P (GeV) 


20 < P T < 70 


13727 


549340 


10 < P T < 140 


22648 


2271090 


Longitudinal P (GeV) 


|-Pl| < 60 


13342 


392401 


\Pl\ < 130 


22459 


2051010 


# of charged tracks 


AT \ in 
Nchd > 1U 


12936 


374877 


AT \ in 
Nchd > 


21270 


1936220 


iviaximum r (^vjev ) 


Prnax ^ 30 


11743 


205038 


Prnax ^ 60 


20556 


1167050 


Y23 value 


Y 23 < 0.02 


7775 


74439 


Y 23 < 0.02 


14992 


465461 


F12 value 


0.2 < y 12 < 0.8 


7438 


62584 


0.2 < Y12 < 0.8 


14500 


413762 


Di-jet mass (GeV) 


100 < Mjj < 130 


6691 


19061 


100 < M j:j < 130 


12334 


71918 


Likelihood ratio 


LR > 0.165 


6293 


10940 


LR > 0.395 


9543 


11092 


Significance (Efficiency) 


S/y/S + B 


47.9 (32.5%) 


S/VS + B 


66.4 (36.3%) 



where Mj 1 j 2 /j. i j 4 , M z / H represent the dijet invariant masses paired from the four jets (j'1-4) and 
the Z and Higgs masses, respectively. Here oz = 4.7 GeV and ou = 4.4 GeV are used for 
y/s = 250 GeV and a z / H = 3.8 GeV for y/s = 350 GeV. They are determined from the dijet 
mass distribution reconstructed from the true MC information. After the jet pairing, background 
reductions are applied. 

To select the four-jet-like events, cuts on the number of charged tracks N c h ar ged and jet clustering 
parameter Y34 are applied. Y34 is the minimum scaled jet mass y required for four-jet clustering. 
The leptonic backgrounds (££££, Uqq) are reduced effectively by these selections. In addition, cuts 
on the thrust and thrust angle are applied to reduce the ZZ background, utilizing the difference 
between the event shape of the signal (spherical) and ZZ, qq (back-to-back). The numbers of qqqq 
and qq background events are reduced by a cut on the angle between the Higgs candidate jets (Oh)- 
The WW and ZZ backgrounds are further suppressed by cuts on the Higgs and Z candidates after 
the kinematical constraint fit is applied to the four-jet system as follows. Each jet is parameterized 
by Ej t , 9{, and fa (i = 1 — 4) and fitted with constraints on the total energy = V*)> the 

total momentum (Yli-^ji = 0), and Higgs and Z mass difference (\Mj 1 j 2 — Mj 3 j 4 \ = \Mh — Mz\), 
where Ej^ Pj v 9i, and fa are the energy, momentum, and theta and phi angles of the i-ih jet, 
respectively. After these cuts are applied, an additional cut is applied on the LR derived from 
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the following input variables: thrust, cos thrust > minimum angle between all jets (6 m in), number of 
particles in Higgs candidate jets, fitted Z mass, and fitted Higgs mass. The likelihood cut position 
is selected to maximize signal significance; LR > 0.375 for 250 GeV and LR > 0.15 for 350 GeV. 
All background reduction procedures are summarized in Table II. The background fractions after 
all cuts are 80% qqqq and 20% qq at 250 GeV and 60% qqqq, 30% qq and 10% it at 350 GeV. 

TABLE II: Summary of qqH channel background reduction assuming C = 250 fb _1 with P(e~,e+) = 
(-0.8, +0.3). 



CM energy (GeV) 


250 


350 


Cut names 


condition 


Sig. 


Bkg. 


condition 


Sig. 


Bkg. 


Generated 




52507 


45904900 




36099 


22210900 


x 2 


x 2 < 10 


32447 


2608980 


x 2 < 10 


20207 


1034810 


# of charged tracks 


N chd > 4 


25281 


1120950 


N chd > 4 


14900 


305649 


I34 value 


- iog(y 34 ) > 2.7 


25065 


1002125 


- io g (y 34 ) > 2.7 


14543 


250995 


thrust 


thrust < 0.9 


24688 


935950 


thrust < 0.85 


13522 


144560 


thrust angle 


|cos0 thrust | < 0.9 


21892 


696201 


COSflthrust < 0.9 


12523 


107025 


Higgs jets angle 


105° < 6 H < 160° 


20062 


622143 


70° <0 H < 120° 


11185 


77659 


Z di-jet mass (GeV) 


80 < M z < 100 


16359 


411863 


80 < M z < 100 


9468 


45671 


H di-jet mass (GeV) 


105 < M H < 130 


16359 


411863 


105 <M H < 130 


9451 


44399 


Likelihood ratio 


LR > 0.375 


13726 


166807 


LR > 0.15 


8686 


25393 


Significance (Efficiency) 


S/^/S + B 


32.3 (26.1%) 


S/VS + B 


47.1 (24.1%) 



C. Leptonic channel (£ + £ H) 

For leptonic channel analysis, we considered the cases where the lepton is an electron or a muon. 
We considered only the ttqq and ivqq background processes. First, the following cuts were applied 
to selected isolated leptons: 

• Lepton isolation: E cone < 20 GeV (cone angle: 10°), 

• Lepton track momentum: 

10 < Eiep < 90 GeV at ^ = 250 GeV, 
10 < Ei^ < 160 GeV at yfs = 350 GeV, 
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where E cone is the energy sum for particles within 10° of the lepton. The prompt lepton has a 
smaller E cone than nonprompt leptons. Electrons and muons are identified from their charged 
tracks as follows: 

. Electron ID: > 0.9, 0.7 < < 1.2 

Erotai P 

• Muon ID: — < 0.5, — - — < 0.4, 

Erotai P 

where Eecal, Erotai and P denote the ECAL energy associated with a track, total energy de- 
posited in the ECAL and HCAL, and track momentum, respectively. If there are more than two 
isolated lepton candidates after the electron or muon identification, a pair whose invariant mass 
is closest to Z is selected. After dilepton identification, forced two-jets clustering is applied to 
the remaining particles and the following selections are applied. First, a dilepton mass (M«) cut, 
which should be consistent with the Z mass, is applied: 70 < Mfj < 110 GeV for electrons and 
70 < Ma < 100 GeV for muons. Because the ZZ or WW backgrounds are boosted to the forward 
region compared to the signal, a cut on the Z direction is applied: |cos#z| < 0.8. Finally, cuts 
on dijet mass (Mjj) and a mass recoil to the lepton pair (M rec ) are applied to select the Higgs 
signal: 100 < Mjj < 140 GeV and 70 < M rec < 140 GeV for electrons; 115 < Mjj < 140 GeV and 
70 < M rec < 140 GeV for muons. The background reduction procedures for the leptonic channel 
are summarized in Table III. After all cuts were applied, the background was dominated by the 
iiqq whereas the vtqq was well suppressed. 

IV. BRANCHING RATIO MEASUREMENT 

After event selection, the measurement accuracies of the Higgs BRs to bb, cc, and gg are evalu- 
ated on the basis of a template fitting to the flavor likeness of the Higgs dijets obtained by using 
the LCFIVertexing package [15]. The probabilities of b and c quarks for each jet q (i = 1,2)] 
are calculated in LCFIVertex using neural net training with a Z — > qq samples at the Z-pole. 
In addition, another c probability (bc\^) is also calculated whose neural-net is trained only with 
Z — > bb sample as the background. For Higgs dijets, we define the flavor likeness X (X = b, c, be) 
as follows from the Xi [xi = bi, Ci, bci (i = 1, 2)] flavor probability of each jet: 

x = X -^2l . ( 2 ) 

X\X2 + (1 - Zl)(l - X 2 ) ' 

The flavor tagging performance in the ZZ — > vvqq sample at the yfs = 250 and 350 GeV is 
shown in Fig. 4. The ZZ — > vvqq samples are compared for each CM energy because they form the 
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TABLE III: Summary of background reduction in the eeH and (i^H channels assuming C — 250 fb 1 with 
P(e~,e + = (-0.8, +0.3)). 



CM energy (GeV) 




250 


350 


Cut names 


e/n 


condition 


Sig. 


Bkg. 


condition 


Sig. 


Bkg. 


Generated 


e 
fj, 




3137 
2917 


4512520 
4512520 




2740 
1789 


3822410 
3822410 


# of e/fi track ID 


e 
M 


iV e >=2 
JV„ >=2 


2717 

2668 


204403 
28175 


n e >= 2 
% >=2 


2270 
1631 


179580 
23598 


Di-lcpton mass (GeV) 


e 
M 


70 < M« < 110 
80 < M u < HO 


2208 
2287 


34162 
12901 


70 < M u < HO 
80 < Mu < 100 


1425 
1406 


51436 
13313 


Z direction 


e 
M 


|cosfl| < 0.8 


1797 
1889 


21600 
8036 


| cos6»| < 0.8 


1192 
1203 


20874 
6250 


'Ol-iot" TTlrtQQ (\Xo\f\ 


M 


100 < Mjj < 140 
115 < Mjj < 140 


1394 
1445 


2721 
1955 


110 < Mjj < 140 
115 < Mjj < 140 


865 
855 


2019 
1197 


Recoil mass (GeV) 


e 
M 


70 < M rec < 140 
70 < M rec < 140 


1184 
1365 


1607 
983 


70 < M rec < 140 
70 < M rec < 140 


567 
638 


590 
465 


Significance (Efficiency) 


e 
M 


S/VS 1 + B 


22.4 (37.8%) 
28.2 (46.8%) 


S/VS + B 


16.7 (20.7%) 
19.2 (35.7%) 



same final state as Z — > qq, which was used to train the flavor tagging neural network. Figure 4 
shows that no significant difference in the flavor tagging performance at y/s = 250 and 350 GeV is 
observed for any of the flavors. 

To evaluate the measurement accuracy of the BRs, the b-, o, and frc-likenesses of the selected 
events were binned in a three-dimensional histogram and fitted with those of the template samples, 
which consist of H — > bb, cc, and gg and other background processes. Figure 5 shows the three- 
dimensional histogram projected to the two-dimensional b- and c-likeness axes for the hadronic 
channel. The probability of entries in each template sample bin is expected to be given by the 
Poisson statistics: 

PiJk = ^ (n = JV$», /i = N^) , (3) 
where and Nf^ a are the probability of entries and the number of data entries at the k) 
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FIG. 4: Flavor tagging performance at CM energies of 250 and 350 GeV in ZZ — > vvqq sample. Horizontal 
axis shows the efficiency of b/c jets; vertical axis shows the purity of tagged b/c jets. 

bin, respectively. jq t ^P late [ 8 given by 

N te r late = £ T s ■ Nf jk + N$ , (4) 
s=bb,cc,gg 

where N?- k is the number of entries at the (i,j,k) bin in each H — > bb, cc, and gg template; 
N^ k 9 is the number of entries in the background template sample, which is the sum of the SM 
background events and the Higgs-to-nonhadronic decay events. Furthermore, r bb , r c5 , and r gg are 
the parameters to be determined by the template fitting. They are defined as the Higgs branching 
fractions to H — > bb, cc and gg, respectively, normalized by that of the SM, 

a-Br(H -> s) 
a SM .Br(H^ s) ' 

Here a is the Higgs production cross section and a SM and Br(H — > s) SM are the cross section and 
branching fraction in the SM, respectively. From Eq. (5), the measurement accuracies of a ■ Br are 
obtained as follows; 



r * = c„ ^ — ZsM ( s = 66 ' c ^ 99)- (5) 



A (a ■ Br) , „ N Ar s , 

'-{H -> s) = — (a = 66, cc, gg). 



a ■ Br y r s 



The r s 's values were determined by a binned log likelihood fitting, where each bin probability is 
given by Eq. (3). 

On the basis of the three-dimensional (3D) histogram, 5000 toy MC events were generated using 
the Poisson distribution function for each bin, which were fitted to obtain r b i, r c5 , and r gg . The 
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FIG. 5: Two-dimensional images of three-dimensional template samples for 6-likeness vs. c-likeness 

number of bins in the 3D histogram were optimized to minimize the statistical fluctuation in the 
fitted results caused by low-statistic bins. Bins with fewer than one entry were not used for the 
fitting. The distributions of r 6 g, r cg , and r gg for template fitting to 1000 toy MC events are shown 
in Fig. 6. The error in r s is determined by the Gaussian fittings to these distributions, which are 
shown in Tables IV and V for CM energies of 250 and 350 GeV, respectively. 

The tables also show the accuracies after correction of the total cross section. From a study of 
the recoil mass in the process of e + e~ — > eeH and HfiH, the accuracy of the total cross section 
(Aa/a) was estimated to be 2.5% at 250 GeV [9, 16]. For 350 GeV, we assumed an accuracy 
of 3.5% because the recoil mass measurement relies on the ZH process, whose cross section is 
inversely proportional to the square of the CM energy; thus, the accuracy of the total cross section 
measurement would be inversely proportional to the CM energy. 

From Tables IV and V, we see that the Higgs cross section times branching ratio can be measured 
at about 1% for H — >• bb and 7 to 9% for H —> cc and gg. The measurement is approximately 
10 - 20% better at 350 GeV than at 250 GeV. The instantaneous luminosity at 350 GeV is 25% 
greater than that at 250 GeV according to the ILC beam parameters. Thus, for an equal running 
time, measurements at 350 GeV will give us about 20 — 30% better accuracy than those at 250 
GeV. On the other hand, the accuracy of the BR to bb, ABr / Br(H — > bb), is limited by the total 
cross section ambiguity; thus, measurement at 250 GeV gives us better results than that at 350 
GeV. In the other decay channels, comparable BR measurements are possible even if the same 
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integrated luminosities are assumed. 




TABLE IV: Summary of template fitting results r s and accuracies of (a ■ Br) and Br after correcting a for 
an accuracy of 2.5% at ^Jl = 250 GeV assuming C = 250 ftr 1 with (e", e+) = (-0.8, +0.3). 





vvB 


qqH 


e+e-H 


H+Ijl-H 


comb. 




1.00±0.02 


1.00±0.01 


1.00±0.04 


1.00±0.03 


1.00±0.01 




1.02±0.11 


1.01±0.10 


1.02±0.27 


1.01±0.23 


1.02±0.07 


T 99 


1.02±0.14 


1.02±0.13 


1.05±0.33 


1.02±0.24 


1.02±0.09 


*^(H^bb) (%) 


1.7 


1.5 


3.8 


3.3 


1.0 




11.2 


10.2 


26.8 


22.6 


6.9 


^l(H^gg) (%) 


13.9 


13.1 


31.3 


33.0 


8.5 


^(H^bb) (%) 


3.0 


2.9 


5.7 


4.5 


2.7 


^(H^cc) (%) 


11.4 


10.5 


31.3 


22.8 


7.3 


^(H^gg) (%) 


14.2 


13.3 


33.1 


24.0 


8.9 



V. CONCLUSION 

The measurement accuracy of the Higgs branching fractions, H — > bb, cc, and gg, were eval- 
uated at y/s = 250 GeV and 350 GeV. In terms of signal significance, yfs = 350 GeV yields 
better background suppression than y/s = 250 GeV for each channel. The combined results for 
measurement accuracies of the Higgs cross section times BRs (A(cr ■ Br) / a ■ Br) to H — > bb, cc, and 
gg are 1.0%, 6.9%, and 8.5% at CM energies of 250 GeV and 1.0%, 6.2%, and 7.3% at 350 GeV, 
respectively, assuming the same integrated luminosity of 250 fb _1 . At the ILC, the total Higgs 
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TABLE V: Summary of template fitting results r 8 and accuracies of (a ■ Br) and Br after correcting a for 
an accuracy of 3.5% at </s = 350 GeV assuming C = 250 fbr 1 with (e~, e+) = (-0.8, +0.3). 





VV £1 


qqti 


c e n 


fl [X n 


comb. 




1.00±0.01 


1.00±0.02 


1.00±0.05 


1.00±0.05 


1.00±0.01 




1.02±0.11 


1.01±0.10 


1.02±0.31 


1.04±0.32 


1.01±0.06 




1.02±0.14 


1.04±0.14 


1.04±0.37 


1.03±0.34 


1.02±0.07 


Ur >(H~>bb) (%) 


1.4 


1.5 


5.3 


5.1 


1.0 


->«•)(%) 


8.6 


10.1 


30.5 


30.9 


6.2 


*<£*p.{H^gg) (%) 


9.2 


13.7 


35.8 


33.0 


7.3 


^L(H^bb) (%) 


3.8 


3.8 


6.4 


6.2 


3.6 


^(F^cc)(%) 


9.2 


10.6 


30.7 


31.1 


7.2 


*^{H^gg) (%) 


9.8 


14.1 


36.0 


33.2 


8.1 



cross-section a is measured using the Z recoil mass process. Using Aa/a = 2.5% for 250 GeV 
and assuming it is 3.5% at 350 GeV, Higgs BRs (ABr/Br) to bb, cc, and gg are derived as 2.7%, 
7.3%, and 8.9% at CM energies of 250 GeV and as 3.6%, 7.2%, and 8.1% at 350 GeV. Therefore, 
we conclude that the Higgs cross section times BR (Br x a) can be measured better at 350 GeV 
than at 250 GeV owing to the higher S/N at the higher energy. However, when the accuracy of 
the total cross section measurement by recoil mass measurement is considered, BR of H — > bb can 
be measured better at 250 GeV, even if the integrated luminosity is the same at both energies. 
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